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Abstract 

A method to study weakly bound three-body quantum systems in two 
dimensions is formulated in coordinate space for short-range potentials. 
Occurrences of spatially extended structures (halos) are investigated. Bor- 
romean systems are shown to exist in two dimensions for a certain class 
of potentials. An extensive numerical investigation shows that a weakly 
bound two-body state gives rise to two weakly bound three-body states, 
a reminiscence of the Efimov effect in three dimensions. The properties 
of these two states in the weak binding limit turn out to be universal. 

PACS number(s): 03.65.Ge, 21.45. +v, 31.15.Ja, 02.60Nm 

Introduction. Characteristic properties of halo systems are binding energies 
much smaller than the typical energy of the interaction and spatial extensions 
much larger than the range of the potential |l], ^ . Three-body halos in three di- 
mensions exhibit many interesting features. Borromean systems, discovered now 
experimentally, are bound three-body structures where none of the two-body 
subsystems are bound [|, The Thomas and Efimov effects are anoma- 
lies due to a singularity which occurs in a three-body system when the ratio 
of scattering length and effective range is infinitely large §, ft |, g Then 
the effective three-body potential behaves as the inverse square of the distance 
which results in infinitely many three-body bound states. 

The Thomas effect corresponds to an infinitely small effective range. The 
Efimov effect is associated with an infinitely large scattering length where the 
wave functions of excited states reside in the tail region of the effective three- 
body potential. Physical examples of these might possibly exist in nature [fill 

The three-body halos in two dimensions are less studied and the present 
investigations |L3| |TJ, [l5[ |l6) are limited in various ways, ft has been proven 
[ p~6| that the number of bound states is always finite. However, the possibility of 
existence of a Borromean state has not been examined and therefore the number 
of bound states is not established. Moreover the properties of these states have 
not been investigated in detail. 
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In this letter we shall for simplicity only consider a system of three identical 
bosons in two dimensions. The purpose is to (i) formulate an efficient method 
to solve the coordinate space Faddeev equations in two dimensions, applicable 
for arbitrary short-range two-body potentials, (ii) derive asymptotic equations 
for the effective three-body potential which alleviate crucially the numerical 
investigations of weakly bound systems, (iii) discuss the possible existence of 
Borromean systems, halos, Efimov and Thomas effects and (iv) illustrate nu- 
merically the universal properties of the three-body halos. 

The method. We shall use the method developed for three dimensions j9j. 
Let Vjk = (vj — Tk) be the distance between particles j and k, r^-M = — r* + 
(rj + rfc)/2 the distance between particle i and the center of mass of particles 
j and k. All particles have the same mass m. The Jacobi coordinates are then 
introduced as Xj = r,-fc/-y/2, y< = r i(j&) -\/2/3. The hyperspherical coordinates 
in two dimensions are given by {p,Qi} = {p,ai,9 xi ,9 y i}, p = \J x\ + yf, an = 
arctan(ccj/yi), 9 X i and 9 y i are the azimuthal angles of Xj and yj. The volume 
element in hyperspherical coordinates is p 3 dpd9 X id9 y i sinc^ cosaidat and the 
kinetic energy operator is 



A 2 = -^-j -2cot(2oi)— (2) 
Oaf oai 

Id 2 Id 2 



sin 2 ( ai )de 2 xl cos 2 (a 4 ) d0 2 yt 



The total wavefunction is now expanded in a complete set of hyperangular 
functions 

1 °° 

*(p,n) = — £/»(p)*»(p,n). (3) 

^ n—1 

where $ n (/0, fi) for each p are chosen as eigenfunctions of the hyperangular part 
of the Schrodinger equation 

(A 2 + ^ J2 VirM $„ = A n (p)$„ . (4) 

Here V is the two-body potential and the expansion coefficients f n (p) satisfy 
the system of coupled equations 

d 2 A n + 3/4 2mE\ 

Q i 2+—^-Qnn-^ r )fn{ P ) (5) 
= fann' +2P nn >^-\ fn-(p) , 
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r 02 
Qnn'(p) = J dQ $* n (p,tt)-^$ n ,(p,n) (6) 

Pnn'(p) = Jdfl **(p,n)^-* n ,(p,n) . (7) 

The wavefunction $„ is written as a sum of three components, each expressed 
in the corresponding system of Jacobi coordinates 

»=i 

these components satisfy the three Faddeev equations 

A 2 0i J) + 2 -^ P 2 V{n)<5> n = A„( P )#, i = 1,2, 3. (9) 
a 

The physical solutions are equivalent to the solutions of the Schrodinger equa- 
tion in Eq.(0) but these equations are better suited for descriptions of subtle 
correlations. 

The s-wave motion is responsible both for the long distance behavior and 
the Thomas and Efimov effects. We therefore restrict to s-waves where the 
wavefunction only depends on p and on. The Faddeev components in Eq.(^|) 
must all be expressed in the same set of Jacobi coordinates. This amounts for s- 
waves to rewriting one set of coordinates in terms of the other and a subsequent 
integration over the angular variables 9 X and 8 y , i.e. 

<f>(p,a') ^ ±- J <f>(p,a'(a,(3))d8, (10) 

where a' is given by a and the variable 8, describing the rotation between two 
sets of Jacobi coordinates, given by 

sin 2 a' — - sin 2 a + - cos 2 a ± sin a cos a cos 8, (11) 
4 4 2 

where the choice of ± is independent of 8. For three identical bosons {4> = 4>n ) 
the Faddeev equations then reduce to the three identical equations 

.^_ 2cot(2a) ^_ A ^ (Pia)= (12) 

2m r- . . / , . , 1 



ft 2 



p 2 V(V2p sin a) (d{p, a) + ^ J 4>{p, a') d(3^j 



3 



Large-distance behavior. We define an angle ao such that 1 2mp 2 V{y/2p sin ao)/h 2 1 
|A(p)| and we assume that the potential is of short range, that is ao approaches 
zero with increasing p. If a > ao the potential is then negligible for large p and 
Eq.(|T2|) becomes 

-^-2cot(2a)^--A^(a)=0. (13) 



da 2 da 

The solution which satisfies the boundary condition of vanishing derivative at 
a = ir/2 is 

4>(a) =7rcos7w P v (cos 2a) — 2 sin -kvQ u (cos 2a), (14) 

where A = Av(v + 1) and P and Q are Legendre functions. At small a this 
solution behaves as 

4>{a) = 2sin(z/7r) (7 + log a + ?/>(l + u)) 

+7rcos(z/7r) + 0{a 2 ) , (15) 

where ip is the digamma function and 7 is Euler's constant. Without interactions 



the solution in Eq.(14) with the boundary condition of zero derivative at a = 
provides the quantization rule v = 0,1, 2..., see Eq.(|l5|). 

For a non-zero short-range potential the integral in Eq. (|l2| ) can be expanded 
for a < a Q -C 1 as 

- 4>{p, a') dp = 2*£) + 0(a 2 ). (16) 
t Jo a 

Then Eq.(Q2h simplifies to a differential equation with an inhomogeneous term 



-\)<l>(a) = -^p 2 V{V2p sin a)2<^(7r/ '3) . (17) 

The homogeneous part of this equation reduces for small a by the substitution 
r = pay/2 to the two dimensional Schrodinger equation for two particles. Thus, 
assuming that A can be neglected compared to 2mp 2 V/fi 2 , the large-distance 
(a < ao <C 1) solutions to the homogeneous part of Eq.(|l7j) simply are the 
zero-energy two-body solutions. Then for a ~ ao the physical solution is ap- 
proximately C log (y/2pa/a) , where C is an arbitrary constant and a is the 
scattering length defined as the distance where the two-body wave function for 
zero energy is zero. One solution to the inhomogeneous part of the equation is 
now — 2</>(7r/3). The complete physical solution to the inhomogeneous equation 
is therefore (a ~ ao) 



( V2pa 



(a) = Clog^^J - 2 ^UJ- ( 18 ) 
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Matching the two solutions, Eqs.(0) and (|l8|), and their derivatives at a = 
ao gives the equation 

2 sin(^7r) log f ^— - ] = 2 sin(f7r) (7 + ip(l + v)) 
\ a ) 

+7T C0${VK) + 2(j) . (19) 

Both for p a and p <C a the logarithm at the left hand side is large. The quan- 
tity v must therefore approach an integer I to compensate for this divergence. 
The leading order of an expansion in powers of p for the lowest v gives 



-1 

3 
2 



log 



V 3a J 



(20) 



Therefore this eigenvalue A = 4^(i/+l) is approaching zero in both these limiting 
cases in contrast to three dimensions where a negative constant asymptotically 
is approached when the scattering length is infinitely large. The effective radial 
potential in Eq.(|5|) has therefore a repulsive centrifugal term and no collapse of 
the wave function in the center is possible. 

Let us now consider diverging solutions \(p) to Eq. (p"9"|) . Then large imag- 
inary values of v oc — ip are necessary and the outer function in Eq.([l4|) ap- 
proaches 

^ )w ^/2£ sin ( ( ^ + ^ )(2a ^ )+ i) ■ (21) 

which for large p is exponentially small at a = ir/3 compared to it's value at 
a = ao. For large imaginary v the equation in Eq.([l9|) then becomes 

=^+^+i-^ +i h (22) 

which has the asymptotic solution 

„• / . /o „ \ 

+ 0( P - 2 ) (23) 
(24) 

This parabolic behavior of A is the signature of a bound two-body state with 
the binding energy B and wave number k given by 4e~ 27 a~ 2 = 2m*B/h 2 = k 2 , 
where m* — m/2 is the reduced mass of the two particles. We can verify this 
by solving the two-body problem, where the s-wave radial solution outside the 
potential is Ko(kr) sa — log (fcr/2) — 7 which must be matched with the solution 
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inside the potential at ro = paov2. For small binding we can use the zero- 
energy solution which at r ~ ro is log (r/a). This matching gives the above B 
and k and the result is accurate to the order ro/a. For attractive zero-range 
potential it is therefore exact. For the potentials without two-body bound state, 
where ro/a is not small, this particular solution does not exist. 

The wavefunction is exponentially small everywhere except in a small region 
close to a — 0, see Eq. (pl|) . In this region the zero-range wave function obtained 
from Eq.(|l7|) is proportional to K (kr = kpaoV2) and after normalization ap- 
proximately given by 

<j)(a) ps 2kpK (V2kpa) . (25) 

The related diagonal part of Q n is then by use of Eq.(||) computed to be Qu — 
~ ip 7 ' wrucn m combination with Eq.(p4]) gives the first diagonal equation in 
Eq.|| as 

/ d 2 1 2m, A „, , 

This is the characteristic large-distance behavior of a two-body radial Schrodinger 
equation in two dimensions. 



Efimov and Thomas effects. The eigenvalue equation for large distances is 
in three dimensions given by § 

where A = v 2 — 4. For p^> a the lowest solution is v — 2 (A = 0) that is a regular 
free solution whereas for p <C a, A — * Aoc = —5.012..., which leads to a strongly 
attractive p~ 2 potential in the radial equation. Physically the p~ 2 behavior 
is limited at small distances by a finite range R of the potential, and at large 
distances by a finite scattering length a. However, in the limit R/a — > such 
a potential gives rise to the so called "falling towards the center" phenomenon 
with infinitely many bound states. When R — > 0, the phenomenon is called the 
Thomas effect and when a — > oo it is called the Efimov effect. 



Eq.(27) is in two dimensions replaced by Eq.(|19|) and the Efimov effect is 
therefore not present in two dimensions. Furthermore, the hyperradial poten- 
tial has even for zero-range potentials a repulsive centrifugal barrier at small 
distance. Then the Thomas collapse is not possible and the three-body system 
must for zero-range potentials have a finite number of bound states. 

Equations (^) and ( |27| ) can be formally written in terms of hypergeometric 
function as one general equation for <i-dimensions, where d is a real number. 
The asymptotic solution Aoo of this equation leads to the Efimov and Thomas 
effects in the region 2.3< d <3.8. 
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Figure 1: The effective radial potential, U — a 2 [(§ + Ai) p)~ 2 — Qui , as function 
of hyperradius (logarithmic scale) for a zero-range two-body potential. The scattering 
length is a. The two unnormalized bound- state wave functions are also plotted. 



Bound states. The solutions to the zero-range potentials are similar to those 
of purely attractive weak potentials. The two-body system has at least one 
bound state and for such interactions Borromean systems do not exist in two 
dimensions. The three-body ground state is more bound and an excited state 
is in addition always present. In Fig.l is shown the effective radial poten- 
tial together with the resulting two bound-state wave functions. The corre- 
sponding root mean square radii are in units of the scattering length respec- 
tively < (r/a) 2 > x / 2 = 0.111,0.927. Their large sizes are reminiscent of the 
three-dimensional Efimov states with an extension comparable to the scattering 
length. 

The energies E3 of these states are given in terms of the two-body bound- 
state energy E 2 as E 3 /E 2 = 16.52, 1.267. We have numerically tried to find a 
third bound state with weaker binding energy and larger radial extension by 
calculating the zero-energy wavefunction and looking at the number of nodes. 
Even a careful search to about 10 3 times the scattering length did not reveal 
another bound state. This strongly indicates that only two bound states exist. 
The large proportionality factor for the ground state energy can be considered 
a reminiscence of the Thomas effect. These relations are still valid for arbitrary 
weakly attractive, finite-range potentials. This is because weak binding corre- 
sponds a scattering length much larger than the range of the potential, which 
is the limit of a zero-range potential. 

In Fig. 2 we show the ratio (E 3 — E 2 )/E 2 as function of \E 2 \ for different 
potentials. We first notice that there are always two states which in the weak 
binding limit approach the results for zero-range potential (marked with small 
circles). The purely attractive potentials (solid lines) as well as repulsive core 
potentials (dashed lines) fall approximately on the same universal curves. This is 
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Figure 2: Ratio of three- to two-body energies as function of the two-body energy 
for different two-body potentials V(r) = [Si exp {-\r 2 /b 2 ) + & exp (-2r 2 /b 2 )] . 

The unspecified strength parameter Si is used to vary the two-body binding. 



equivalent to the observation in three dimensions that only low-energy scattering 
properties like the scattering length are important in the weak binding limit. 

The potentials with a short-range repulsive barrier (dash-dot and dotted 
lines), unlike the repulsive core, produce energies deviating in the middle of the 
plot from the universal curves. The reason is that a sufficiently large barrier 
and a sufficiently large attraction produce a three-body ground state confined 
inside the barrier, that yields three bound states in total. As j^l then de- 
creases towards zero two cases are possible. If the potential does not allow the 
spatially confined ground state the third bound state disappears and the first 
two approach the zero-range limit (dotted lines). If, however, the potential al- 
lows a spatially confined ground state, all three states survive in the small \E 2 \ 
limit (dash-dot lines). Now the second and the third bound states approach 
the zero-range limit, while the ground state energy remains finite. In this case 
the ground state persists even into the region where two-body state is unbound, 
creating thus a Borromcan phenomenon. This is only possible with a repulsive 
confining barrier which therefore also limits the spatial extension of the three- 
body system. The properties of Borromean systems are therefore different in 
two and three dimensions. 

Conclusions. Based on the hyperspherical expansion of the Faddeev equa- 
tions we have investigated the possible structure of three weakly bound iden- 
tical bosons in two dimensions. For purely attractive two-body potentials and 
the potentials with repulsive cores the Borromean systems do not exist. For 
two-body potentials with a short-range repulsive barrier the two-body system 
may not have a bound state while a three-body bound state exists. Borromean 
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systems are therefore possible in two dimensions. 

We find numerically that a weakly bound two-body state is always accom- 
panied by two bound three-body states, resembling the Efimov states in three 
dimension. If a Borromean state is present there are therefore in total three 
bound three-body states. For all types of potentials, two of these three-body 
bound states have energies and radii following a universal curve in the weak 
binding limit. Their sizes scale with the two-body scattering length and can 
therefore become arbitrarily large in analogy to s-state halos in three dimen- 
sions. 
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